arXiv:1502.01960vl [math.PR] 6 Feb 2015 


Collective periodicity in mean-field models of cooperative behavior 

F. Collet'^, P. Dai Pra^, M. Formentin'^ 


Dipartimento di Matematica, Alma Mater Studiorum Universitd di Bologna, Piazza di Porta San 

Donato 5, 1-40126 Bologna, Italy 

^Dipartimento di Matematica, Universitd degli Studi di Padova, Via Trieste 63, 1-35121 Padova, Italy 
‘^UTIA, Czech Academy of Sciences, Pod Voddrenskou vezi 4, CZ-18208 Prague, Czech Republic 


Abstract 

We propose a way to break symmetry in stochastic dynamics by introducing a dissipation 
term. We show in a specific mean-field model, that if the reversible model undergoes a phase 
transition of ferromagnetic type, then its dissipative counterpart exhibits periodic orbits in 
the thermodynamic limit. 


1. Introduction 

Many real systems comprised by many interacting components, for example neuronal net¬ 
works, may exhibit collective periodic behavior even though single components have no 
“natural” tendency of behaving periodically. Various stylized models have been proposed 
to capture the essence of this phenomenon; we refer the reader to [8] and [6j, where many 
related references are given, and possible origins of this periodic behavior are discussed in 
clear and rigorous terms. Most of the proposed models, despite of their rather simple struc¬ 
ture, turn out to be hard to analyze in rigorous terms; their study ends up to the search of 
stable attractors of nonlinear, infinite dimensional dynamical systems. 

The main purpose of this paper is to illustrate and study a class of models that are obtained 
as perturbations of “classical” reversible models by adding a dissipation term, which breaks 
reversibility. Models of this sort have already appeared in different context (see e.g. m 
for applications to multicellular dynamics, and [7] for applications to finance). The simplest 
interacting system in this class is the dissipative Curie-Weiss model proposed and analyzed 
in pp. Its infinite volume limit can be reduced to a two dimensional nonlinear dynamics, 
whose large-time behavior can be determined explicitly. It is shown, in particular, that by 
tuning the parameters of the model a globally stable periodic orbit emerges through a Hopf 
bifurcation. Note that this phenomenon is qualitatively similar to the emergence of periodic 
orbit in a system of interacting FitzHugh-Nagumo models presented in [8]. 

In this paper, after having presented the general procedure to add dissipation in a system, 
we concentrate on the dissipative version of a model for cooperative behavior which, in its 
reversible version, has been extensively studied in [2] and, more recently, in [5]. This model, 
characterized by a quartic self-potential and a mean-field interaction, shares many features 
of interacting FitzHugh-Nagumo models. We show that the thermodynamic limit can be 
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studied in a suitable small-noise approximation, which reveals various peculiar features that 
we find interesting. 

• A stable periodic orbit arises through a homoclinic bifurcation rather than a Hopf 
bifurcation; this is a global phenomenon, that is not detected by a local analysis. 

• The stable periodic orbit may coexist with stable fixed points: the long-time behavior 
depends on the initial condition. 

• The system can be excited by noise: by increasing the noise, stable fixed points may be 
de-stabilized, and trajectories are attracted by the periodic orbit which can be quite 
far from the fixed point. 

We remark that the small noise approximations we obtain are quite similar to the class of 
Gaussian nonlinear processes originally studied in m and [12], and more recently in m- 
In particular, im and |12] give an interesting overview of oscillatory phenomena in these 
systems, and of the role of the noise. 

In Section we describe a general class of dissipative models we consider. In Section we 
derive the macroscopic limit of a specific class of models for cooperative behavior. The be¬ 
havior of the macroscopic dynamics is studied in Section]^ in particular, we study in details 
the system in absence of noise and then we consider a suitable small noise approximation 
of the macroscopic equation, for which we prove noise-excitability. Most proofs are then 
postponed to Section and to the Appendix. 

2. Systems driven by a random potential 
2.1. Reference model 

Let x(f) = denote the positions at time t oi N particles moving in Particles are 

subject to a potential field, consisting of two parts: the external potential and the interaction 
potential. The external potential is a smooth function [/ : —)■ M. Moreover, particles 

themselves generate a potential field, which is felt by all other particles. We adopt the 
mean-field assumption that particles contribute symmetrically to the field. More specifically, 
letting be the space of Borel probabilities of provided with the weak topology, 

and denoting by 

1 ^ 

P7v(x(f)) := — e 

i=l 

the empirical measure of the particles at time t, we assume the potential at x G generated 
by the particles to be of the form 


Vt{x) = V(x, pjv(x(t))) (2.1) 

for some sufficiently regular function P : x —)> M. Finally, for a given family 

{wi{t) : t > 0)^1 of independent Brownian motions, we consider the evolution given by the 
following system of stochastic differential equations 

dxi{t) =—VU{xi{t))dt— Wt{xi{t))dt + adwi{f) (for f = 1,..., A"). (2.2) 
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2.2. Introducing dissipation and diffusion 


In the model above, the interaction potential in (2.1) is a deterministic function of particle 
positions. When dissipative effects are present, as in a great variety of models in biology 
(see e.g. [3 H [131 IS]) this scheme should be modihed, by letting the potential Vt evolving 
according to the following (Stochastic) Partial Differential Equation 


dVt{x) = —aVt{x)dt + DAVt{x) + dV{x, /9Ar(x(t))) + noise, (2-3) 


where a,D > 0 tune dissipation and diffusion, respectively. The noise term, which could 
have several different forms, will be removed in the examples below. 

Remark 1. A similar construction can be applied to systems with discrete state space, e.g. 
spin systems for which Xi G { — 1,1}. In this case the Laplacian is omitted in the evolution 
) of V), and the evolution \2.2\ is replaced by a Glauber dynamics for the time-dependent 
potential Yhi {U{xi) + Vt{xi)) at some inverse temperature /3. The special model correspond¬ 
ing to U = 0 and 

- I ^ 

V(x, pn(x)) = -X— ^ Xi 

i=l 

is studied in [1]. 



3. A model of cooperative behavior 


We consider here d = 1 and the following choice of external and interaction potential: 

4 2 


and 


0 


N 


V{x,pn{x)) =—2_^{x - Xi) , 
i=l 

where 0 is a positive parameter that represents the strength of the interaction between 
particles. The corresponding reference (non-dissipative) model has been extensively studied 
in [2] . More recently, the same model has been applied to describe systemic risk in finance 
(see [5]). It is known that in the limit as N ^ -|-oo the system exhibits a phase transition: for 
small 6 the mean particle position is zero, while for large 9 particles tend to cluster around 
the two minima of the double-well potential U{x), giving rise to “polarized” equilibria. 


To avoid complications, we assume the initial condition Vo{x) to equation (2.3) is a polyno¬ 
mial in x: 


n 


Vq{x) := > afc(0)x 


fe=o 


Observing that, by Ito’s rule 


dV (x, /97v(x(t))) = 9xd 



N 

i=l 


N 


N 2 

'^Xi{t) dxi(t) + —dt, 


i=l 
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equation (2.3) is solved by the polynomial 


Vt{x) := '^ak{t)x^ 


k=0 


whose coefficients satisfy the system of equations 
dak{t) = —aak{t)dt 

< daklt) = -aaklt)dt + D{k + 2){k + l)ak+ 2 {t)dt 
dai{t) = —aai{t)dt + 6Da3{t)dt + 9d 


for k = n,n — 1 

for 2 < k < n — 2 (3 1) 


Note that ao(t) is not needed for (2.2). By (3.1), it follows that ak{t) —)■ 0 as t —)• +oo, for 
all k > 2. Since we are interested in a steady state regime, we may assume V)(x) is linear 
in X for all t. Renaming fi{t) := ai(f), we obtain the following form for ( |2.2| ) coupled with 

(pl): 


dxi{t) = (—+ Xi{t)) dt — fj,{t)dt + adwi{t) (for i = 1,..., N) 

dfi{t) = —afi{t)dt — 6dm^^\t), 

where we set := ^ Note that ( |3.2[ ) can be rewritten as 

dxi{t) = (^—x^{t) + Xi{t)) dt — fj,{t)dt + adwi{t) (for i = 1,..., iV) 

dd{t) = -{(^-d)fk{t)dt-+ 


(3.2) 


(3.3) 


We remark that, although the drift terms in (3.3) are not uniformly Lipschitz, strong ex¬ 
istence and uniqueness can be established, for example by using the classical Hasminskii’s 
Test: by defining, for instance 


1 


N 


nx,/i) 

2 = 1 


4 2 ' 

4 2 


Qj 9 

+ 2 ^ 


for a > 0 sufficiently small, one obtains an inequality of the form 

CV{x,fi) < K[l + V{x,fi)] 


for some K > 0, where £ is the infinitesimal generator of the diffusion (3.3). This inequality 
implies existence and uniqueness of strong solutions (see e.g. M)- 

3.1. Propagation of chaos 

In the limit as N ^ +oo, the system of equations (3.3) naturally suggest the following 
limiting process, which describes the behavior of a single component in the limit: 

dx{t) = (—x^(t) -|- x{t)'j dt — fi{t)dt + adB{t) 


= -(a - 9)fi(t) + BE [x^(t) - x(t)] , 


(3.4) 


where B{-) is a standard Brownian motion. Note that this equation has a non-/oca/structure, 
due to the appearance of the law of x{t) in the term E [x^(t) — x{t)] of the right hand side. 
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Theorem 2. For each /io G K and each real random variable with finite third moment and 
independent of the Brownian motion B{-), equation (3.4) with initial condition x(0) = 
/i(0) = fj,o, has a unique strong solution. 


Theorem 3. Let be the solution of 

isfying the following conditions: 

(a) are independent and identically distributed with a given law A having finite 

third moment; moreover, they are independent of the Brownian motions {wi{-))^i. 


(3.3) with an initial condition sat- 


(h) ^W(O) = /xo G M. 

Then, as N ^ +oo and for each fixed k > 1, the vector-valued process 




dm 


im. 


restricted to a given time interval [0,T] converges in law to (yi(-)) y2(')) • • ■ > where 

yi(')) 2/2(')) • • •) yA:(') independent copies of the solution of ( |3.4| ) with initial condition 
yj{0) ~ A, ^(0) = ^ 0 - 

The proofs of Theoremsandare rather standard; for completeness, we sketch them in the 
Appendix. Note that the only difficulty is to deal with the non-global Lipschitz property of 
the drift. In the microscopic model this problem is overcome by stopping the process at the 
boundary of a compact set and controlling this stopping time via Lyapunov methods. This 
approach is not directly applicable at the macroscopic level: stopping the process affects the 
drift globally, due to the non-locality of the evolution. 


4. Main results 


Our main goal is to study the long-time evolution of solutions to equation (3.4) and of the 
corresponding Fokker-Planck equation: 


dtqtix) = Y^xxQtix) - dx {[-x^ + X - nit)] qtix)} 
^(t) = -{a-e)p.{t)-9{-x^ + x,qt) , 


where the notation (/, q) 
derivative guarantees that. 


:= f f{x)q{x)dx is used. The regularizing effect of the second 
for f > 0, the law of x{t) has a density qt{x) which solves (4.1). 


f.l. Stationary solutions 

The study of equilibria for the macroscopic dynamics turns out to be essentially trivial. 


Proposition 4. Equation (4.1) admits a unique stationary solution given by 

.,.4 


q^{x) = • exp ^ 


1 




-Y + ^ 


/X* = 0 


(4.2) 


where is a normalizing factor. 


We remark that the scenario is quite different than for the reference model with the same 
potential but without dissipation, where multiple equilibria arise for large 9. 
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J^.2. The noiseless dynamics 

Letting cr = 0 in (3.4), we obtain a deterministic dynamics described by the ODE 


_(t) = + x{t) - n{t) 

= e)ix{t) + e [x^{t) - x{t)] . 


(4.3) 


The dependence of the attractors for (4.3) on the parameters is quite nontrivial, despite of 
the fact that (4.3) has always (0,0), (—1,0) and (1,0) as unique equilibrium points in the 
plane (x,/i). Moreover, for all values of the parameters a > 0 and 0 > 0, the equilibrium 
(0,0) is a saddle point. We denote by W the stable manifold of the origin, i.e. the set of 
initial conditions whose corresponding solution converges to the origin as t —>■ +oo. 


Theorem 5. Fix a > 0. Then there is 9i = 0i(a) > 0 such that the following properties 
hold. 

(a) For 6 < 9i, W is an unbounded, open curve, that separates the basins of attraction of 
(—1,0) and (1,0) respectively. 


(b) For 9 = 9i, W closes into an eight-shaped curve, comprised by two cycles, symmet¬ 
ric with respect to the origin, surrounding (—1,0) and (1,0) respeetively (homoclinic 
bifurcation). 


(c) For 9i < 9 < a -\- 2, W spirals around two unstable periodic solutions of (4.3), sur¬ 
rounding (—1,0) and (1,0) respectively. These cycles separate the basins of attraction 
of (—1,0) and (1,0) from the basin of attraction of a unique stable periodic solution 


of (4.3), which surrounds both (—1,0) and (1,0). 


(d) At 9 = a 2, the two unstable cycles eollapse in (—1,0) and (1,0) through a Hopf 


bifurcation. For 9 > a + 2, (—1,0) and (1,0) become unstable: (4.3) admits a unique 
periodic solution, whose basin of attraetion is \ [{(±1,0)} U W]. 


4 . 3 . Excitability by noise: Gaussian approximation 


The study of the effects of the noise {a > 0) in (3.4) appears to be of a higher level of difficulty 


compared to the analysis of the deterministic system. Simulations of the interacting particle 


system dynamics (3.2), indicate that the following picture arises. 


The structure of the attractors is preserved for small a, except that the fixed points 


(±1,0) are replaced by the invariant law (4.2). For a sufficiently large, the noise breaks 
the organized structure, so that, in particular, there is no periodic solution. 

For intermediate values of a the behavior depends on the parameters. The most in¬ 
teresting case is for 9i < 9 < a ± 2. By increasing a, the basin of attraction of the 


invariant law (4.2) shrinks, and this invariant law becomes unstable at a critical value 

> (Tc. Thus, an 


(T = (Tc. A unique stable periodic solution survives for a < a(, with 
initial condition that, for small a, would be in the basin of attraction of (4.2), gets 
exeited by the noise, and is attracted by a periodic orbit. See Figure [TJ 
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Figure 1 : The figures illustrate the phenomenon of excitation by noise, that is how noise (of suitable 
intensity) helps/favors the emergence of a stable periodic behavior. Simulations of the interacting 
particle system (|3.2|) have been run for 9 slightly below the threshold value a + 2. Recall that this 


range of the parameter corresponds to a coexistence phase for the deterministic system (4.3), in 
which the locally stable equilibria (±1, 0) coexist with a macroscopic stable limit cycle. We simulate 
the trajectory t i— (t), fj,{t)) for different noise sizes. Both simulations are performed for the 
same number of iterations, with the same time-step size, and starting from /i(0)) = (0,0) 

with particles initially placed half in (-1-1,0) and half in (—1,0). Only the noise size is changing. On 
the left panel, when a is small, the curve {m^^\t), fi{t)) stays close to the initial point for all times; 
indeed, it oscillates with maximum amplitude of order of 10“®. On the right panel, we increase the 
noise size and (m^^^(t),/i(t)) exhibits a periodic behavior. 


We have no proof of the above statements, but strong evidences is obtained through numer¬ 
ical simulations. Surprisingly, even the analysis of the linear stability of (4.2) appears to be 
quite hard. One possible approach consists in considering the infinite dimensional system of 
equations for the moments of (3.4). By Ito’s rule we obtain, for k > 1: 

^ r . .. . ^2 

dx^{t) = 


—kx^~^‘^{t) + kx^{t) + ^k{k — l)x^ ‘^{t) — kfj,{t)x^ ^{t) 


dt + akx^ ^{t)dB{t). 


Setting mk{t) := E(x^(t)), taking the expectation in the previous equation we get 
druk cr'^ 

—^{t) = -kmk+ 2 {t) + krukit) + - l)mk- 2 {t) - k^i{t)mk-i{t) 

^(t) = -{a - 9)n{t) - 6 + mi{t)) = -aix{t) - 


(4.4) 


This system of equations appears to be as hard as (4.1). However, we use it as a guide to 
construct a Gaussian approximation of (3.4). To be precise, denote by {xa{t), t^a{t)) the 
solution of (3.4), for a given deterministic initial condition X(j(0) = x,;U(j(0) = /U (random 
initial conditions could be considered as well). We construct, on the same probability space, 
a process {yo-{t),i 2 o-{t)) satisfying the following conditions: 

(a) ya{-) is a Gauss-Markov process, Vai') is a deterministic process, ycr(O) = x, zZo-(O) = y. 


(b) The moments of ya{-) and zZo-(-) satisfy (4.4) for A: = 1,2. 

(c) For every T > 0 there is a constant Ct such that for all cr > 0 

E 


sup \Xait) - yait)\ 
te[o,r] 


-b sup \ya{t) - w(A)| < 
te[o,r] 


(4.5) 
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Note that (4.5), which expresses the fact that is a small noise approximation of Xo-, is not 
enough to identify the law of y^, since this condition in insensitive to corrections of order cj^. 
These corrections are determined by condition (b); indeed, the equations (4.4) for k = 1,2 
form a closed system of two equations, since for Gaussian random variables, all moments 
can be expressed in terms of the first two. 

We now define the Gaussian process y^- First, consider the following system of ODEs: 


drrin 


= -mi + ma - Va - 3a^maVa 


dt 

dVfj dm„ 

■dT = 

^^ = l + 2{l-3mJ)C,-6<72Vj 

ma{0) = X = k- K(0) = 0. 


(4.6) 


It is easily seen that solutions of (4.6) are bounded, which implies global existence and 
uniqueness. We can therefore define the following Gauss-Markov, centered process, as the 
unique solution of the linear SDE: 


dzait) = [l — — 3cT^V'o-(t)] Zcr{t)dt + dB{t) 

2 :^( 0 ) = 0 , 


(4.7) 


where B{t) is the same Brownian motion driving (3.4). By computing one checks that 

Var{za{t)) = Va{t). 


Finally, let 


ya{t) := mail) + aZa{t). 


(4.8) 


Proposition 6. Properties (a), (b) and (c) hold for the process y^ defined in (4.8) 


Having obtained a small noise Gaussian approximation for (3.4), we study the approximated 


process. Being Gaussian, it is enough to study the evolution of mean and variance. Since 
E(yo-(t)) = ma{t) and Var{ya^{t)) = cj^I 4 -(t), we are left to the analysis of the three dimen¬ 
sional system (|4.6|). In the following analysis we omit the subscript a in mo-, t'o-, W-- In the 


three dimensional space {(m, iZjV) : m, o € M., V > 0}, equations (4.6) admit the following 
equilibrium points: 


f + Vl — 3cr^ ^ 1 — Vl — 3it^ 

•si = I -- 


V2 


6 cr 2 


f ^/l + Vl — 3(7^ ^ 1 — Vl — 3cJ^ 

S2 = I-- 


V2 


6 cr 2 


( Vl - Vl - 3^ 1 -h Vl -3 (t2' 

•53 = I-7^^-, 0 , - 


V 2 

for V < 1/3, and 


6(7^ 


S5 = 0, 0, 


sa = 


' Vl - Vw^^ 1 + Vl - 3cj2' 


V2 


6 cr 2 


1 -|- Vl + 6 (t2 

6 cj 2 


for all CT > 0. Note that, as cr 0, the equilibria si and S 2 converge to (±1, 0, 3/2), so they 


correspond to the equilibria (±1,0) of the deterministic system (4.3); the equilibria S 3 , S 4 
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and S 5 converge to (0, 0, +oo). In the following result we give a local version of the property 
of excitability by noise: under suitable conditions on the parameters, by increasing a the 
equilibria si and S 2 lose stability. 


Theorem 7. Assume 9 < a+ 2, with a+ 2 —9 sufficiently small, and assume a < 10. Then 
there exists Uc > 0 such that si and S 2 are linearly stable for a < Uc, but linear stability is 
lost as a crosses ac- 


For 9 > a + 2 and u <C 1, we expect system (4.6) to have an unique attracting periodic 
solution, as the deterministic system (4.3) has. This behavior is illustrated by simulations 
in Figure 




- 2-1012 t 


Figure 2: Regime cr <C 1 and 0 > a + 2. In the left panel the (Ft = 0)-section of the typical motion 
of (toct, Ucr) is shown. The trajectory 1 1 — (jna-ft), Va{t)) (solid red line) approaches and remains close 
to the limit cycle attractor for the dynamics (4.31 (dashed black line). A heuristic explanation of the 
phenomenon can be inferred from the right panel, where the oscillatory behavior of me and Ft is 
illustrated. Since Ft is of order one for all times, the terms proportional to in (4.6) are small 
when cr ^ 1. Thus, the trajectories of (4.6) are close to those of the deterministic system (4.3). 


The behavior of (4.6) for large a could be studied as well, although it may not reflect any of 
the features of (3.4). For 9 > a + 2 one shows, in particular, that the equilibrium s^ becomes 


linearly stable as cr > y |(q; — 0)(q! — 0 — 1), but it is not a global attractor, since a locally 
stable periodic orbit survives. This is illustrated in Figure 


5. Proofs 

5.1 . Proof of Proposition 

In view of dynamics B, an equilibrium probability density for our system must satisfy 
the system of equations 


cr 


- 9a: { + x - /i*] g*(x)} = 0 


— {a — 9)fi^ — 9 (—+ X, g*) = 0 , 


(5.1a) 

(5.1b) 


Solving (5.1a) we obtain 


q^{x) = ^ • exp <j ^ 


—^—h X + 20/x*x 
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Figure 3: Simulations for the system (4.61 in the case when 9 > a+2 and a > ^— 9){a — 9 — 1). 
In this regime the limiting behavior of the trajectories depends on the initial value of the magne¬ 
tization. We run simulations with starting point Fct(O) = Va-ifi) — 0 and various initial conditions 
for TOct(O). For |mo-(0)| <C 1 the oscillatory behavior disappears. Indeed, (m^, F^.) converges to 

S 5 = (0,0, ), which is the only fixed point for the dynamics (4.6) in this phase and is linearly 

stable (left panel). On the contrary, when 1771^(0)1 is large, the periodicity persists also for cr 1 
(right panel). 


with 



1 




—^—h -|- 29^^x 


dx. 


This result together with equation (5.1b) imply 


0 and the proof is complete. 


5.2. Proof of Theorem\^ 

The scenario depicted in the statement is due to the occurrence of both a homoclinic and a 
Hopf bifurcation. To ease the readability of the proof, we believe it is worth first to explain 
what is going on and then to give the technical details. 

Roughly speaking, the theorem states there exist three possible phases for system (4.3): 

• Fixed points phase. For 9 < 9i the only stable attractors are (±1,0). 


• Coexistence phase. For 9 = 9i the system has an eight-shaped, symmetric with respect 
to the origin, homoclinic orbit surrounding (±1,0). By increasing the parameter 9 
from 9i, this separatrix splits into an outer stable limit cycle, which contains both 
the equilibria, and two inner unstable periodic orbits around ( 1 , 0 ) and (— 1 , 0 ), re¬ 
spectively. In this phase (±1,0) are linearly stable. Therefore, two locally stable fixed 
points coexist with a stable limit cycle. 


• Periodic orbit phase. For 9 = a + 2 a Hopf bifurcation occurs: the inner unstable 
limit cycles disappear collapsing at (±1,0). At the same time the equilibria lose their 
stability and, thus, the external stable limit cycle remains the only stable attractor 
for 0 > a ± 2 . 


The key point of the proof is the particular structure of the vector field generated by (4.3). 
For convenience, let us denote such a vector field by 


Va,9{x,ia) 


( —x^ + x — fj, 

y—(a — 9)pL ± 9{x^ 



(5.2) 
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Observe that (5.2) defines a one-parameter family of negatively rotated vector fields (with 
respect to 9, for fixed a); that is, the following are satisfied 


the critical points of (5.2) are isolated 


at ordinary points, as the parameter 9 increases, all the field vectors Vap{x,iJ,) rotate 
clockwise. 


Additionally, (5.2) is semi-complete] in other words, the field rotates of an angle vr as the 


parameter 9 varies over M. 

For dynamical systems that depend on this specific way on a parameter, many results con¬ 
cerning bifurcations, stability and global behavior of limit cycles and separatrix curves are 
known. For the sake of completeness and readability, we collect here the properties satisfied 


by (5.2) that are crucial in the sequel, referring to m Chapter 4] for precise and general 
statements. 

(1) Limit cycles expand or contract monotonically as the parameter 9 varies in a fixed 
sense. 

(2) A limit cycle is generated/absorbed either by a critical point or by a separatrix of 


(3) Cycles of distinct fields do not intersect. 


Now let us go into the details of the proof. 

Lyapunov function. Let a,9 > 0. Consider the function 

{9x + p.y 


W{x,^^) = - + 


2a9 


and observe that the total derivative 


^W{x,pL) = -x'^ + (1 + 9)x^ - \{9x + /i)^ 

at u 

is negative for every n and sufficiently large x. Thus, there exists a stable domain for the 
flux of (4.3) and, in particular, the trajectories can not escape to infinity as t —)• +00. 


Local analysis of equilibria and Hopf bifurcation. As already mentioned, the vector field (5.2) 


admits three fixed points in the phase plane (x, p,): (0, 0) and (±1, 0). The origin is a saddle 
point for all values of the parameters; while, the equilibria (±1,0) are linearly stable for 
0 < a ± 2 and their local stability is lost for 9 > a ± 2. At the critical point 0 = a ± 2 a 
Hopf bifurcation occurs: by decreasing the parameter 9 from a ± 2 two symmetric unstable 
limit cycles are generated from (±1,0). This proves the first assertion in statement (d). 

Separatrix formation and stability. Properties [T] and allow us to explain the appearance 
of a separatrix curve, whose breakdown causes a homoclinic bifurcation at 9 = 9i. Indeed, 
while decreasing 9 from a ± 2, the cycles arisen from (±1,0) expand until they join each 
other at the origin forming a homoclinic eight-shaped separatrix graphic Tq at 9 = 9i. Notice 
that, for 9 = 9i, Tq = W. This proves statement (b). 

Furthermore the exterior of the separatrix is stable. In fact, the existence of a stable domain 


for the flux of (4.3) implies that all trajectories in an outer neighborhood of Fq approach Fq 
itself as t —>■ ±00. 
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Description of the phases. We are left to prove statements (a), (c) and part of (d). Hence, 


(c) When the separatrix graphic Fq splits increasing 9 from 9i, it generates a periodic 
orbit surrounding both (±1,0) on its exterior and two smaller limit cycles, around 
(—1,0) and (1,0) respectively, on its interior. 

The inner cycles are unstable (due to the subcritical Hopf bifurcation at 0 = a ± 2) 
and represent the boundaries of the basins of attraction of (±1,0). Moreover, the 
external periodic orbit inherits the stability of the exterior of Fq and so it is stable. 
See Theorem 2 in m Section 4.6] for more details. 


(a) It suffices to prove that in this phase the dynamical system (4.3) does not admit a 
periodic solution. Indeed, the non-existence of cycles together with the existence of a 
stable domain for the flux of (4.3) guarantee that every trajectory must converge to an 
equilibrium as t —?• ±00. The fact that the basins of attraction of (±1,0) are separated 
by VF is a consequence of the Stable Manifold theorem (see lini Section 2.7]). 

Thus, it remains to show that there is no limit cycle for 9 < 9i. From propertiesand 
[^it follows that, as 9 increases from 9i to infinity, the outer stable limit cycle expands 
and its motion covers the whole region external to Fq. Similarly, the two unstable 
cycles contract from the graphic Fq and terminate at the critical points (±1,0). As a 
consequence, for 9 > 9i the entire phase space is covered by expanding or contracting 
limit cycles. Now, by using property]^ we can deduce that no periodic trajectory may 
exist for 0 < 0i. In fact, such an orbit would intersect some of the cycles present when 
0 > 01 that is not possible. 


(d) The statement easily follows by combining the impossibility of escaping trajectories, 
the presence of the stable manifold of the origin, the loss of stability of the fixed points 
(±1,0) and the existence of a stable limit cycle as a unique attractor in the phase space. 


To conclude, we observe that 0i > 0. If 0i was negative, for 0 = 0 some periodic solution 
should exist. It is sufficient to show that this is not the case. For 0 = 0, the trajectory of /r 
in (4.3) is /i(f) = /r(0) and clearly excludes the possibility of having limit cycles. 


5. 3. Proof of Proposition 

Let Ztj{t) be the solution of the linear stochastic differential equation (4.7). Clearly Z(j{-) is a 
centered Gaussian process, so ya{t) '■= mfj{t) + aza{t) is obviously Gaussian, which establish 
property (a). Moreover E(yo-(f)) = m^^^t), and Var{yfj{t)) = a‘^Vcr{t). Property (b), i.e. the 


fact that the hrst two moments of ya{t) and t'o-(-) satisfy (4.4) for k = 1,2, follows from a 


simple direct computation, that uses equations (4.6), and the fact that 


IE [yl{t)] = ml{t) ± 3a‘^m„{t)V^{t), E [yt{t)] = m^(t) ± 6a‘^ml{t)V„{t) ± 3a^V^{t). 


So, we are left with the proof of Property (c). By the first equation in (4.6) and (4.7), we 
have 


[ Q o o n 

—m); ± — 3a rUaVa — 3amfza ± az^ — 3a Va-Za-j dt ± adB 

= [-yl + ya + 3a‘^ma{zl - I 4 ) ± a'^{zl - 314^^^)] dt ± adB. 


(5.3) 
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By subtracting ( |5.3[ ) from (3.4), we get 

Xa{t) - ya{t) = [ { [yl{s) - xl{s) - y^is) + Xo-(s)] + K(s) “ y^a{s)]} ds -cr^ [ Ra{s)ds 
Jo Jo 


= [Xa{s)-ya{s)][l-f{s)]ds+ [uais)-lJ.ais)]ds-a'^ Ra{s)ds, 

Jo Jo Jo 


where 

and 


Rait) := 3m^(t)(4(t) - T4(t)) + cr{zl{t) - 3Va{t)z„{t)) 


fid) ■= xlit) + yl{t) + x„{t)y„{t) > 0 . 


(5.4) 

(5.5) 


Note that (5.4) can be seen as a linear differential equation for := Xa{t) — yait) of the 
form ^ ^ 

'!/;(t) = [ A{s)'tlj{s)ds + J B{s)ds, 


/o 


with initial condition V’(O) = 0, whose solution is given by 


ifit) = / i3(s)exp 

Jo 


A{r)dr 


ds. 


(5.6) 


This gives 


Xa{t) - yait) = f e f‘>d /(^))«*»’ [u^{s) - iJais) - a'^Rais)] ds. (5.7) 

Jo 

On the other hand, by subtracting from the equation ^i^ait) = —avait) — 0^E(yo-(t)) the 
corresponding equation for /i^, we obtain 

T^ait) - yait) = -a [ [i^ais) - yais)]ds - 9[E{y^{t)) - E(Xcr(t))]. 

Jo 

This last formula provides an equation for (t) — ya it) of the form 

V’(t) = [ A{s)'ip{s)ds + C{t), 


whose solution is 


ip{t) = C{t)+ / ^(5)0(5) exp 

Jo 


A{r)dr 


ds. 


This yields 

I'ait) - yait) = 6 
from which the following inequality follows: 


(5.8) 


E{xait)) -E{y^{t)) - a / e *)[E(x^(s)) - E(y^(s))]ds 

Jo 


sup \vait) - yait)\ < 0 sup |E(Xf^(t)) - E(y(^(t))| 
te[o,T] te[o,r] 

— aO sup f |E(xo-(s)) — E(yo-('S))| ds 

te[o,r] Jo 

<Kt sup E[\xait) - yait)\], 
i6[0,T] 


(5.9) 
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for some constant Kt > 0. By (5.7) and (5.9) we get 


E 

sup \xa{t) - ya{t)\ 

< At [ E 

sup \x„{s) - ya{s)\ 


i6[0,T] 

Jo 

se[o, 4 ] 


ds + Ato-"^ sup E[|i?o-( 0 |] 
te[o,r] 

for some At > 0 which, by Gronwall’s Lemma, yields, for all cr in a bounded subset of 
( 0 , +oo) 


E 


sup \Xa{t) - ya{t)\ 
i6[0,T] 


< Btct^ sup E[|i?o-(t)|] < Ct(7^, 

te[o,T] 


(5.10) 


where the last inequality follows from the fact that Ra{t) is a polynomial function of a 
Gauss-Markov process and, therefore, has a L^-norm which is locally bounded in time. To 
complete the proof of (4.5) we still need to show that 

sup lucrit) - Z/cr(t)| < 
te[o,r] 


for some Ct > 0, which follows immediately from (5.10) and (5.9). 

5.4- Proof of Theorem^ 

By symmetry, it is enough to consider the equilibrium si. Gonsider the vector field V^q 
associated to equation (4.6) 

( —m? + m — v — Sa'^mv 
— {a — 6)v + 9{m^ — m) + 39a‘^mv 
1 + 2(1 — 2,m?)v — 

We linearize the system around si; we obtain the linear system associated to the Jacobian 
matrix 


/ 


(^i) = 


-1 - Vl - 3 ( t 2 -1 


3 a 


1 + Vl — 3(T^ ) —a + 9 —2>9a‘^ 
3\/2 


2\/rTvr^3^ \ 

+ Vl — 3(^ 


^/2 


V \/l + \/l — 3cj2 


In particular, for cj^ = 0 and 9 = a + 2, 


0 


-3- 

-2 

-1 

0 

2 (a + 2 ) 

2 

0 

3 

0 

-4 


whose spectrum is given by 

Spec (si)] = {-4} U {±iV^} ■ 

Denote by A(cj) the continuous function giving an eigenvalue of the matrix DV ^,^_^_2 {si), 
with A(0) = \f2ai. We are going to show that, if a < 10, then 


(i(cj^ 


-Re (A((t)) 


> 0 . 


(7 = 0 


(5.11) 
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It follows that there is cr > 0 such that the matrix DV^^j ^2 ('^i) has an eigenvalue with 
strictly positive real part. By continuity, this implies that if 0 < a + 2 but sufficiently close 
to a + 2 , then also the matrix DV^q (si) has an eigenvalue with strictly positive real part, 
and this suffices to complete the proof. Thus, we are left with the proof of (|5.11 ). We write 


0 = det ('Si) - A(fT)l) = ao(cr) + ai(fT)A((T) + 02 ( 0 -) (A((t))^ + 03 ( 0 -) (A(fT))^ , 

which gives 


d(cr^ 


(A(ct)) 


d 


(T=o ai(0)d((T2) 


ao(fT) 


(T=0 


Since 


det (I?F ,%+2 (si) - AI) 


= -A'^ - 2 


1 + Vl - 3u2j + [2 - a + 12 cj 2 - (a + 2) Vl - 3^2 


A 


we obtain 


d 


d{a‘^ 


(A(a)) 


cr=0 


Inserting A(0) = we easily get 

d 


d{a‘^ 


Re (A((t)) 


— Aa\/\ — ^1 + \/l — 3(7^j , 

3A2(0)+ 3 (5 +f) A(0) + 18a 
3A2(0) + 8A(0) + 2a ‘ 

3(10-a) 


cr=o 2(8 + cr) 


from which (5.11) follows. 


6. Appendix 

6.1. Proof of Theorem^ 

It will be convenient to write {|3.4|) in the form 


dx{t) = (—x^(t) + x(t)) dt — fj.{t)dt + adB{t) 
^(t) = - 0^E [x(t)]. 


( 6 . 1 ) 


The second equation in (6.1) can be solved in 

/i(t) = e“"Vo — 0IE[x(t)] + 0e“"*E(^) + a9 f (g_2) 

Jo 


To show existence and uniqueness in (6.1) we follow the argument in [2] (Theorem 2.4.1). 


Given /cq £ 1^; we define the stochastic processes {xnit))t>o and the functions {Hn{t))t>o, for 
> 0, by the following Picard iteration: 

dxn{t) = (-x^(t) + Xn{t))dt - Hn{t)dt + adB{t) 

Hn+iit) = - 6E[xn{t)]+6e~°‘^E{^) + aO [ e~°‘J~''''E[xn{s)]ds, ^ ^ 

Jo 
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all with the same initial condition Xn(0) = /in(0) = /io- Setting 

fit) := xl^^it) + xlit) + Xn+lit)Xnit) > 0, 

and observing that 

xl+i{t) - xlit) = ixn+lit) - Xnit))fit), 


we get 


Xn+lit) - Xnit) = [Xn+lis) - Xnis)][l - fis)]ds - [(J-n+lis) -/J-nis)] ds 

Jo Jo 

= [Xn+iis) - Xnis)][l - fis)]ds + 9 E (s) - x„_i(s)] ds 

Jo Jo 


t ps 


— a9 



e '’^E [x„(r) — x„_i(r)] drds. 


0 JO 


Note that, setting V'(i) := Xn+iit) — Xnit), as in the proof of Proposition]^ this last identity 
is of the form of a linear differential equation 


ipit) = f Ais)il^is)ds + f Bis)ds, 

Jo Jo 


whose solution is given by (5.6). This yields 

Xn+lit) - Xnit) 


= / e->» 

Jo 

Now, set 




6E[xnis) - Xn-iis)] - a9 [ e '’^E [xn(r) - x„_i(r)] dr 

Jo 


ds. (6.4) 


ipnit):= sup E [|Xn+i(s) - X„(s) 
se[o,t] 


It follows readily from (6.4) that for each T > 0 there is a constant Ct > 0 for which 


^n+lit) <Ct (pnis)ds, 

Jo 


which implies 


^niT) < ifliT) 


iCrT) 


n—1 


(n — 1)! 

As a consequence, the sequence of functions E[xn(-)] is Cauchy in C([0,T]), and therefore it 
converges to a limit m(t). Let x(-) be the unique solution of 

dx(t) = (—x^(t) + x(t)) dt — fiit)dt + adBit) 

pt fa c::\ 

/i(t) = e““Vo — Omit) + de““^E(^) + a6 

Jo 

By mimicking the above argument one gets 


Xn+l(t) - x(t) 

Jo 


dE [x„(s) — x(s)] — ad / e ^^E [x„(r) — x(r)] dr 


ds, (6.6) 
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where, now, 


fit) := ^n+l (i) + X^(t) + Xn+lit)x{t) > 0. 

As before, this implies that E(xn(t)) —>• E(x(t)), and thus m{t) = K{x{t)). This proves that 


x(-) is a solution of (6.1). The very same argument is used to prove uniqueness. Indeed, if 
y{-) is another solution, we write the integral equation for y{t) — x{t) and show as before 


that E{x{t)) = E{y{t)) for all t. Thus y(-) and x{-) are both solutions of (6.5) with the same 
m(-) and the same initial conditions, from which x(-) = y(-) follows. 

6.2. Proof of Theorem\^ 

We follow a coupling method; we refer to [H] for more details on this approach to prove 
propagation of chaos. Let i^i)fLi be independent random variables with law A, and indepen¬ 
dent of the Brownian motions (tCj(f) : t > 0)fLi- Moreover, let {xi^\-), be the 


solution of (3.3) with an initial condition x) TO) = y) TO) = /xoi while (yj(-)) ^(•))i=i 
solve (3.4) with the same initial conditions and B{-) = Wi{-). The claim in Theorem is 
proved if we show that (the apex N in x\^'^ is omitted in what follows) 


lim E 

A^^+oo 


sup \xi{t) - yi{t)\ 
te[o,T] 


= 0 . 


(6.7) 


The strategy is analogous to that of the proof of Theorem We first write, for f{t) := 
Xiit) + yfit) + xi{t)yi{t) > 0: 


xi{t)-yiit)= [ [xi(s) - 7/i(s)][l -/(s)](is - f [y{s) - u{s)]ds, 
Jo Jo 


yielding 


xi{t)-yi{t)= / [y{s) - u{s)]exp 

Jo 


[1 - fir)]dr 


ds, 


which gives, for all t G [0,T], 


\xi{t)-yi{t)\<CT lyis) - u{s)\ds 

Jo 


for some constant Ct > 0. In particular, we have 


E 


sup |xi(t) - yi(t)| 
te[o,r] 


< Ct / E[|t('S) — i^(s)|]ds. 

Jo 


On the other hand, it holds 


( 6 . 8 ) 


(6.9) 


/i(t) — vit) = —a / [/i('S) — r'is)]ds 
Jo 

N 

I 

- e 


1 ^ / 1 ^ \ 
— ^x*(t) - E[yi{t)] - — ^Xi(O) -E[yi(0)] 

' i=l i=l J 


( 6 . 10 ) 
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Set 


N / N 

nYI “ IE[yi(t)] - f ^ X] “ IE[yi( 0 )] 


i=l 


i=l 


By (5.8), 


n{t) — ^{t) = —9ip{t) +a6 ^{s)e "*-* ^^ds. 


( 6 . 11 ) 


Note that if in the expression of ip{t) we add and subtract ^ ~ y*(0)]) the following 

estimate follows for s G [0,T]; 

N N 

E[bis)!! < ^- 9‘WII + 


N ^ u yu ^ 

2=1 2=1 


+ E 


1 ^ 1 r 1 ^ 

j;^'^yiis)-E[yi{s)] +E —'^yi{0)-E[yi 


2=1 


2=1 


< E 


sup \xi{r) - yi{r)\ 

re[0,s] 


+ 


c 

7^’ 


( 6 . 12 ) 


for some C > 0, where the last inequality follows from the fact that E[|xj(t) — yi{t)\] does 
not depend on i, and that the random variables yi{s) are i.i.d. and with finite variance. By 
(6.12) and (6.11) it follows that 


E[|/i(t) - iy{t)\] < At 



sup |xi(t) 
te[o,r] 


yi{t)\ 



(6.13) 


for some strictly positive constant At- Inserting (6.13) in (6.9) we obtain, for some possibly 
different constant Ct and for all t G [0,T]: 


E 

sup Xl(s) 

-yi(s)l 

< Ct / E 

sup \xi{r) - yi{r)\ 


se[o,t] 


Jo 

rS[0,s] 


ds 


Ct 

77’ 


which implies 


E 


sup |xi(t) 
i6[0,T] 


yiit)\ 


< 


Kt 

Tiv 


for some Kt > 0, from which the conclusion follows. 


(6.14) 
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